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ABSTRACT 

Mean  humidity,  temperature,  wind  data  and  fluctuating 
temperature  and  velocities  data  obtained  from  the  R/V  Acania 

were  interpreted  on  the  basis  of  empirically  derived  expres- 

.  .  2 

sions  describing  the  temperature  structure  parameter,  CT  , 

and  the  dissipation  rate  of  turbulent  kinetic  energy,  e  . 

2  .  . 

Properties  of  CT   and  z    are  related  to  a  stability  parameter 

Z/L  and  the  Richardson  number,  Ri .   Non-dimensional  tempera- 

9    l±  I  Q   3  fl   9 

ture  structure  functions  parameters,  Cp  /Z    (-5— )  ,  and  tur- 

-3 
bulent  kinetic  energy  dissipation  rates,  Ks  ZU...   ,  versus  Ri 

were  examined  for  open  ocean  conditions. 

This  experiment  verified  earlier  hypotheses  that  the  pro- 
files obtained  from  empirically  derived  expressions  using 
overwater  data  are  very  much  different  from  those  obtained 
using  overland  data.   The  distribution  of  the  non-dimensional 
turbulent  kinetic  energy  dissipation  rate  agreed  with  that 
predicted  by  similarity  theory. 
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I.   INTRODUCTION 

The  marine  boundary  layer  is  that  region  in  the  first 
kilometer  over  the  sea  which  is  defined  principally  by  tur- 
bulent transfer  of  momentum  and  heat.   In  the  past  decade 
there  have  been  concentrated  efforts  to  make  measurements 
in  this  region  and  to  define  physical  expressions  describing 
the  turbulent  intensities.   The  marine  environment  presents 
a  number  of  problems,  especially  in  accessibility,  platform 
stability,  and  sensor  protection.   Most  important  of  these 
problems  are  the  lack  of  an  extensive  data  bases  and  the 
fact  that  available  turbulence  expressions,  which  are  empiri- 
cal in  nature,  have  been  formulated  primarily  from  overland 
data. 

With  the  development  of  systems  dependent  on  laser  tech- 
nology, the  Department  of  Defense  interest  in  turbulence  in 
this  region  pertains  to  assessing  and  predicting  the  environ- 
mental influence  on  optical  wave  propagation.   The  refractive 
index  problem  for  the  marine  boundary  layer  is  more  compli- 
cated than  the  similar  problem  over  land  because  of  the  dif- 
ficulties in  making  point  measurements  of  the  fluctuating  air 
temperature  from  which  refractive  index  statistics  for  optical 
waves  are  often  obtained.   For  optical  waves  the  atmospheric 
index  of  refraction  variations  primary  result  from  tempera- 
ture-induced density  changes.   These  variations  cause  scintil- 
lation and  phase  fluctuations  of  laser  beams  propagated  in  the 
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atmosphere  and  the  loss  of  sight  of  objects  viewed  through 
high  magnification  telescopes. 

An  interdisciplinary  group  at  the  Naval  Postgraduate 
School  has  been  involved  in  overwater  experiments  to  relate 
the  optical  properties  and  the  optically  relevant  meteoro- 
logical parameters  of  the  marine  boundary  layer  to  observ- 
able bulk  atmospheric  parameters. 

The  primary  purpose  of  this  work  was  to  evaluate  pre- 
viously formulated  overland  expressions  for  the  overwater 

regime.   Expressions  presented  by  Wyngaard  et  al.  (1971)  and 

2 
Friehe  (1976)  were  used  to  relate  C^      to  a  stability  para- 
meter, either  the  Richardson  number  or  Z/L  where  L  is  Monin- 
Obukhov  Length.   In  addition  to  this,  expressions  described 
by  Busch  and  Panofsky  (1968)  were  used  to  relate  e  to  the 
stability  parameters. 

Ten  minute  time  averaged  values  of  wind,  temperature, 
and  relative  humidity  measured  at  four  levels  (4.2,  6.6, 
11.3,  and  17.7  meters)  were  examined  in  this  study.   The 
most  recent  and  the  majority  of  the  observations  examined 
were  made  during  the  months  of  September  and  October  19  76  off 
the  coast  of  San  Diego,  California.   Data  previously  evalua- 
ted by  Hughes  (1976)  and  Karach  (1976)  were  incorporated  into 
the  final  evaluations  to  supplement  the  data  base.   The  latter 
data  were  also  measured  at  four  levels  (4.2,  6.6,  7.6,  and 
13 . 9  meters ) . 
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II.   THEORETICAL  BACKGROUND 

A.   GENERAL 

2 

CT   and  e  are  optically  relevant  turbulence  parameters  be- 
cause they  define  the  intensity  of  fluctuations  in  the  refrac- 
tive index  and  the  size  or  scale  of  refractive  index  inhomo- 
geneities ,  respectively.   The  scales  of  the  temperature 
inhomogeneities  that  distort  image  resolution  and  degrade 
laser  beam  propagation  in  the  atmospheric  surface  layer  lie 
in  the  inertial  subrange  where  fluctuations  are  isotropic. 

The  latter  characteristic  allows  theories  by  Tatarski 
(1961)  describing  electromagnetic  propagation  to  be  applied 
in  analyses.   Tatarski  considered  the  refractive-index-struc- 
ture parameter, 

r  2  _  [n(x)  -  n(x+r)] 2  .    . 

n  271  KL) 

r 

which  is  related  to  the  temperature  structure  function  para- 
meter,  CT  ,  as 

2  -fi  p   -?    2 

CTZ  =  (79  x  10  b  t-)    l    C  Z  (2) 

T  T^      n 


where  P  is  the  pressure  and  T  the  temperature  of  the  atmos- 

2 
phere.   CT   may  be  evaluated  in  a  number  of  methods. 

The  dissipation  rate  s  defines  the  micro-scale,  I    ,  as 


,3  x. 

O       '  £ 


l_       =    (^  )%  (3) 
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where  y    is  the  kinematic  molecular  viscosity.   The  micro- 
scale  appears  in  optical  parameters  for  image  resolution. 
For  example,  the  exp(-2)  folding  distance  of  the  mutual 
coherence,  p   is  defined  as 

p   =  1 (4) 

°    (1.5C  I    "1/6  K) 
n   o 

where  K  is  the  average  wave  number. 

B.   TEMPERATURE  STRUCTURE-FUNCTION  PARAMETER,  C^ 

The  temperature  structure-function  parameter  may  be 
evaluated  or  estimated  by  several  different  methods  depend- 
ing on  the  measurement  procedures.   At  separations  r  of  the 
order  of  inertial  subrange  scales,  the  temperature  structure- 
function  in  a  locally  isotropic  field  has  the  form 

r    2  -  <  [TOO  -  T(x+r)]2>  ,    * 

LT  271  K     J 

r 

2 
Cm   can  also  be  related  to  the  one-dimensional  tempera- 
ture spectrum  (The  Fourier  transform  of  the  correlation  func- 
tion with  separations  in  the  streamwise  direction)  in  the 
inertial  subrange  as 

ST(K1)  =  .25  CT2  K~5/3  (6) 


K   is  the  streamwise  component  of  wavenumber. 


16 


A  third  approach  is  based  on  the  hypothesis  of  Corrsin 
(19  51)  that  in  the  inertial  subrange  the  one  dimensional  tem- 
perature spectrum  can  be  expressed  as 

ST(K)  =  ax£~1/3  K"5/3  (7) 

where  x  is  the  rate  of  molecular  destruction  of  <(T-T)>   ,_ 

av/2 

(in  °C  S"  ),  s  is  the  rate  of  molecular  dissipation  of  turbu- 

2-3 
lent  kinetic  energy  (in  cm   s   ) ,  and  a    is  hypothesized  to  be 

a  universal,  dimensionless  constant;  the  same  in  all  flows. 

2  . 
Equations  (6)  and  (7)  yield  an  expression  for  CT   m  terms  of 

the  rates  of  dissipation  of  turbulent  kinetic  energy  and 

temperature  variance,  viz. 

CT2  =  3xs"1/3  (3) 


where  8  is  an  empirical  constant  with  a  value  of  3.20. 

C.   MICROSCALE  I       (OR  e) 

o 

Both  s  and  x  in  Equation  (8)  represent  molecular  processes 
and  forces  which  act  only  over  very  small  distances.   As  pre- 
viously indicated,  the  particular  scale  length  of  interest 
in  propagation  analyses  is  the  microscale  of  turbulence  de- 
fined by  Equation  (6)  and  is  dependent  on  only  one  dynamic 
parameter,  the  dissipation  rate  of  turbulent  kinetic  energy, 
e,  which  can  be  estimated  from  the  one  dimensional  spectral 
density  of  velocity  as 

-2/3   -S/3 
SU(K)  =  C  £      K  5/J  (9) 

where  E,    is  an  empirical  constant  determined  to  be  0.51. 
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D.  TAYLOR'S  HYPOTHESIS 

Because  measurements  are  normally  made  at  a  point  a  rela- 
tion between  temporal  (f)  and  wavenumber  or  spatial  frequency 
(k)  is  required.   According  to  Lin  (1953),  if  one  considers 
that  eddies  of  wavenumber  K  are  correlated  only  over  distances 
not  large  compared  to  27i7k,  then  over  such  distances  the  non- 
uniformity  of  convection  velocity  should  be  small  relative 
to  the  convection  velocities.   Further,  if  the  mean  velocity 
is  a  good  approximation  to  the  convection  velocity,  then 
0  >>  U1 where  0'  is  the  spatial  gradient  of  U.   If  these  cri- 
teria are  met  the  turbulence  field  may  be  considered  as  being 
frozen  and  swept  at  the  mean  wind  speed  U  past  a  stationary 
instrument.   This  is  known  as  Taylor's  hypothesis  and  it 

allows  frequency  spectra  to  be  interpreted  as  streamwise 

2 

wavenumber  spectra.   This  enables  the  evaluation  of  CT   from 

Equation  (6)  using  point  temperature  fluctuation  measurements. 

E.  RELATING  C^2   AND  E  TO  STABILITY  PARAMETERS 

Turbulence  similarity  theory  relates  both  s  and  x  to 

•  •  2 

measurable  stability  parameters  so  that  CT   can  be  estimated 

using  Equation  (8).   The  1968  Air  Force  Cambridge  Research 
Laboratory  (AFCRL)  boundary  layer  study  yielded  information 
on  the  behavior  of  both  e  and  x  and  the  combination  of  dimen- 
sional and  physical  arguments  that  make  up  surface-layer 
similarity  theory.   The  AFCRL  results  demonstrated  that  the 
statistics  of  both  the  mean  and  turbulent  fields,  when  pro- 
perly non-dimensionalized ,  are  universal  functions  of  the 
stability  parameter  (Z/L). 
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Non-dimensionalized  e  and  x  values  can  be  expressed  as 


— q  =  f,  (Z/L)  (10) 


=  fn  (Z/L)  (11) 


2       2 


where  the  scaling  velocity  is  U...  =  (-uV)  ,  the  scaling  tern- 

perature  is  T.,.  =  ^— n —  ,  the  Monin  Obukhov  stability  length 

is  L  =  -  T  U.,3/kg  wr1-  and  f,(Z/L)  and  f«(Z/L)  are  empirical- 
o  **  1  l 

ly  determined  functions. 

.  .  .         .      2 

An  empirical  expression  relating  CT   to  boundary  fluxes 

arises  from  direct  substitution  of  Equations  (10)  and  (11) 
into  Equation  (8)  as 

2      2   -2/3 
CTZ  =  TV  z      f3  (Z/L)  (12) 


where  f-(Z/L)  is  a  combination  of  f,(Z/L)  and  f„(Z/L).   The 
formulation  of  f-(Z/L)  was  examined  by  Wyngaard,  et  al.  (1971) 
and  their  results  are  shown  in  Figure  1.   They  obtained  the 
following  forms  for  f-(Z/L) 

f3(Z/L)  =  4.9  [l-7(Z/L)]"2/c   0  >  Z/L     (13) 


f3(Z/L)    =    4.9  [1+2.75CZ/D]  0    <    Z/L  (14) 


where  negative  and  positive  Z/L  values  correspond  to  unstable 
and  stable  conditions  respectively. 
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Fig.  1.   The  vertical  profile  of  CT' 


An  empirical  expression  relating  e  to  boundary  fluxes  is 

Z     Z 
Equation  (10),  where  4>m(t-)  ~  f"  can  replace  f,(Z/L),  viz. 


eZK 

u,3 


=  (j>M(Z/L)  -  Z/L 


(15) 


where  <t>M(Z/L)  is  the  non-dimensional  wind  shear.   Overwater 

results  by  Garrett  (1972)  on  the  non-dimensional  dissipation 

_3 
rate  5  KeZU...   ,  versus  Z/L  appear  in  Figure  2. 

Both  Equations  (12)  and  (15)  can  be  rewritten  so  that 

2 
gradient  terms  appear  in  the  non-dimensional  CT   and  z    terms 

if  the  following  expressions  for  T.,.  and  U...  are  used 


KZ     dQ_ 
<J>h(Z/L)  dz 


KZ 


8u 
4>MCZ/L)  Tz 


(16) 
(17) 
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Fig.  2.   The  dimensionless  dissipation  rate  versus  Z/L. 
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Substituting  Equation  (16)  into  Equation  (12)  and  Equation 

(17)  into  Equation  (15)  yield  such  expressions  for  the  non- 

2 

dimensional  CT   and  e  terms  respectively 


2 
C  2 

f„(Z/L)  (18) 


v4/3,86.2      2rZ^    3 
Z    (^}     *h  (L} 


K2[^A(h-2-    I  <tu(h-3]  (19) 


2,8u.3       iyMvL'     L  rMvL 


F.   RELATING  Z/L  TO  RICHARDSON  NUMBER 

The  above  expressions  can  be  interpreted  more  readily  by 

relating  Z/L  to  the  Richardson  number, 

,89. 
Ri  =    .  2  , 


as 


T(H> 


Z/L  =  f^CRi)  (20) 


where  fr(Ri)  is  an  empirical  function  which  has  been  well 
documented  for  the  constant  flux  layer  overland  (Businger 
et  al.,  1971)  and  0   is  the  vertical  potential  temperature 
9v  =  6(1  +  .61q) . 

The  following  relationships  between  Z/L  and  Ri  were 
determined  by  Dyer  and  Hicks  (1970)  for  the  unstable  case 
and  by  Webb  (1970)  for  the  stable  case: 

Z/L    =    Ri  Z/L    <    0 

Ri  (21) 

Z/L    =    -,       P.  Z/L    >    0 

1-aRi 
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Here  a  is  an  empirically  derived  constant  equal  to  5.   These 
relationships  have  been  recently  verified  for  the  overwater 
regime  by  Paulson,  et  al.  (1972). 

Combining  Equation  (18)  and  Equation  (20)  yields  an  ex- 

2  . 
pression  for  CT   in  terms  of  gradient  measurements, 


"47T7TeT2  "  fe(Rl) 


(22) 


The  function  fr(Ri)  which  has  been  determined  for  the  overland 
b 

constant  flux  layer  is  shown  in  Figure  3  (Wyngaard,  et  al., 
1971). 
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Fig.  3.   The  dimensionless  temperature-structure  parameter 
parameter  versus  the  Richardson  number. 
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Similarly,  combining  Equation  (21)  and  Equation  (19) 
yields  the  following  gradient  relationships  for  e 


72,8u.3 


=  K  [(1-15RD2  -Ri(l-15Ri)  '  ]   Ri  <  0         (23) 


=  K2[{l+4.7  C-Si-r)}  2 


72,Su.3      LL^  "       vl-5Ri 


-(1^rI^1+4-7(I?IrI)}"3] 


Ri  >  0         (24) 

The  curve  defined  by  these  relationships  will  be  presented 
later  (Figure  22). 

G.   BULK  AERODYNAMIC  FORMULA 

An  alternate  approach  for  indirect  estimates  of  the  tur- 
bulent fluxes  in  overwater  observational  experiments  is  the 
use  of  bulk  aerodynamic  formulae.   A  bulk  aerodynamic  formula 
is  one  which  relates  the  boundary  fluxes  to  the  wind  speed  at 
a  level  and  the  temperature  and  humidity  differences  between 
that  level  and  the  surface.   Such  formulae  require  several 
assumptions  regarding  the  stability  conditions  of  the  boundary 
layer  and  the  turbulent  processes  within  it. 
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2 

Bulk  aerodynamic  formulae  for  momentum  (U.,.  )  ,  sensible 


heat  (wTf )  and  humidity  (wq ' )  fluxes  have  the  following 
forms , 


U:,2    =    CD   D2  (25) 

w'8    =    C„    OCT    -T    )  =    C„    DAT                         (26) 

n            as  n 

wV   =    Cr    0(Q    -Qo)  =    C^    UAQ                         (27) 

ht           as  Li 


where  the  symbols  have  been  defined  on  page  11. 

Friehe  (1976)  summarized  considerable  amounts  of  recent 
overwater  data  to  obtain  overwater  estimates  of  C~ ,  C„,  and 
Cp.   He  observed  Cp.  and  Cp  to  be  dependent  on  the  values  of 
UAT  as  well  as  height  z.   Friehe  used  the  bulk  aerodynamic 
formulae  to  formulate  an  alternate  expression  for  the  estima- 
tion of  Z/L  in  the  atmospheric  boundary  layer  above  the  sea. 
The  formulation  consisted  of  replacing  the  transfer  terms 
(T...  and  U..c)  in  the  expression  for  L  by  bulk  parameters 
(U.,.2  =  CDU1Q2  and  T.,.)  to  obtain 


f-   =    -,231[^5  +    0.212   ^S]         -20    <    UAT    <    25(mKs"1) 


(28) 
j-   =    -.371[^I  +    0.132    5M]  OAT    >    25(mKs~1) 

L  uz  uz 

Equation  (28)  leads  to  a  new  stability  parameter  X,  defined 

aS  X  =  ^SCI+O.212  ^)   -20  <  UAT  <  25(mKs_1) 

0 

(  29) 
X'    =    MI(l+0.132   $)       OAT    >    25(mKs_1) 
0 
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Note  that  the  signs  of  X  and  X'  are  the  same  as 
Z/L. 

Friehe  defined  X  and  X'  for  the  10  meter  level.   However, 
these  can  be  generalized  to  any  level  by  the  approximation 

Cn(z)     C„(z)     Cr(z)      "VZ  ;  0 

D      -   H      _   E      _   o_  2      r   , 

cn(io)  "  c„(io)  "  cTTIoT   LT  ,z   .J     K6UJ 

D         H         E         LN(Z~} 

o 

This  was  done  in  this  study. 

Substituting  Equations  (2  8)  and  (29)  into  Equation  (12) 

yields  for  -20  <  UAT  <  25(mKs~1) 

CT2  Z2/3      3.12  x  10"3(1+1.62X)"2/3  unstable 


(AT)2        3.12  x  10'3(1-0.635X)    stable 


for  UAT  >  25(mKs~1) 


cT2  Z2/3 


(AT)2 


-3  -2/3       (31) 

.03  x  10  J(1+2.60X')  Z/J      K6±J 
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III.       EXPERIMENT 

A.  PLATFORM  AND  LOCATION 

Observations  aboard  the  R/V  Acania  during  the  experiment 
in  September  and  October  19  7  6  were  made  in  the  area  from  3  2 
15fN  to  34°13'N  and  123°18'W  to  120°20'W;  Figure  4.   This 
area  provided  a  good  sample  for  open  ocean  condition. 

B.  INSTRUMENTATION 

Ten  minute  averages  of  temperature,  wind,  and  relative 
humidity  were  taken  at  four  levels,  on  two  separate  masts 
spatially  separated  on  the  forward  deck  of  the  ship,  Figure  5 

1.  Recording 

Data  logging  during  the  experiments  was  accomplished 
using  the  NPS  developed  MIDAS  (Microprogrammable  Integrated 
Data  Acquisition  System).   This  micro-processor  based  data 
acquisition  system  utilizes  an  Intel  8008  central  processor 
to  control  the  sampling,  averaging,  and  recording  of  mean 
meteorological  data.   All  software  is  written  in  PL/M  to 
facilitate  the  writing  of  self-documenting  program. 

2 .  Mean  Measurement  Systems 

Mean  wind  measurements  were  made  with  a  Thornthwaite 
Associates  cup  anemometer  wind  profile  register  system,  model 
number  104,  Figure  6.   In  operation,  the  shaft  of  a  three  cup 
anemometer  unit  serves  as  the  shutter  between  a  light  source 
and  a  photocell  for  each  revolution.   These  three  plastic 
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Fig.  5.   Mounting  arrangements  aboard  the  Acania, 
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Figo  6.   Thornthwaite  Associates  Cup  anemometer 
wind  profile  register  system. 
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cone  cups  are  reinforced  with  aluminum  frames  and  are  mounted 
at  120°  intervals  on  stainless  steel  tubes  attached  to  the 
main  shaft.   A  characteristic  of  these  cup  anemometers  is  low 
starting  speeds  with  a  small  amount  of  internal  friction  which 
aids  in  checking  inertial  overshoot. 

Mean  temperature  measurements  were  made  with  tempera- 
ture sensitive  quartz  crystal  probes,  (Hewlett  Packard  model 
HP-2850)  at  each  level.   (Figure  7)  RF  signals  from  the 
crystal  probes  and  from  a  reference  oscillator  were  mixed  in 
the  HP-2801A  readout  unit  to  produce  a  beat  frequency  whose 
signature  can  be  analyzed  to  within  0.001  degrees  centigrade 
per  hertz.   Each  sensor  simultaneously  received  pre-experiment 
calibration  against  a  platinum  resistance  wire  thermometer 
in  a  temperature  controlled  circulating  water  bath  over  the 
expected  temperature  range.   The  accuracy  in  achieving  a  0.005 
degree  centigrade  correction  factor  was  a  constant  for  each 
probe.   A  3.7  meter  flexible  coaxial  cable  is  permanently 
attached  to  the  sensor  head  and  the  mast  mounted  probes  are 
housed  in  an  aspirated  shelter  as  seen  in  Figure  8.   Tempera- 
ture values  were  automatically  recorded  on  a  printed  tape. 

Mean  relative  humidity  was  obtained  using  the  Dunmore- 
lithium  chloride  sensor  and  Hygrodynamics  Digital  I  hygrometer 
indicator.   (Figure  9)  The  sensor  was  also  placed  in  the 
aspirated  shelter.   The  basis  of  operations  of  this  system  is 
resistance  charge  within  an  electrolytic  solution  generating 
a  reference  voltage  variance  proportional  to  the  relative 
humidity  change.   Automatic  temperature  compensation  in  the 
instruments  meet  the  following  specifications 
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Fig.  7.   Hewlett  Packard  model  HP- 2 850  temperature 
sensitive  quartz  crystal  probes. 


Fig.  8.   Aspirated  shelter 


Fig.  9.   Dunmore  type 

lithium  chloride 
humidity  sensor. 
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+_  3%  relative  humidity  below  90% 
+_  ■+%  relative  humidity  below  90% 
Calibration  of  sensors  was  accomplished  by  a  comparative 
method  using  a  saturated  saline  solution  in  a  closed  con- 
tainer.  Observation  and  recording  was  accomplished  in  like 
manner  as  the  temperature. 

3 .   Fluctuation  Measurements 

The  temperature  fluctuations  were  measured  using  a 
bridge  developed  at  GTE  Sylvania,  the  GTE  Sylvania  Model  140. 
It  was  necessary  to  modify  the  system  slightly  for  use 
aboard  the  R/V  Acania. 

This  system  is  a  small  temperature  difference  sensing 
device  originally  designed  to  be  carried  aloft  by  a  balloon 
and  equipped  to  process  and  transmit,  by  pulse-rate  modula- 
tion on  a  radio  frequency  carrier,  this  information  to  a 
ground-base  receiving  station  for  demodulation  and  recording. 

The  baseband  portion  of  this  system  is  a  balanced 
Wheatstone  bridge  excited  by  a  3KHZ  signal  with  a  synchronous 
detector  on  the  output.   Segments  of  a  very  small  diameter 
platinum  wire  serve  as  temperature  sensors  in  opposite  arms 
of  the  bridge.   In  the  single  wire  mode,  one  wire  is  replaced 
by  fixed  resistance.   An  output  from  the  bridge  which  is 
proportional  to  the  temperature  difference  between  the  two 
probes  or  a  temperature  change  of  a  single  wire  in  the  single 

wire  mode  is  produced. 

—  fi 

The  sensor  wire  is  0.5  cm  long  and  2.5  x  10   m  in 

diameter.   This  extremely  small  mass  allows  a  response  to 
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temperature  variations  of  up  to  1  KHZ  while  electronic  ampli- 
fication allows  temperature  difference  as  small  as  0.004 
degrees  C  to  be  observed. 

Wind  fluctuations  were  recorded  on  a  Honeywell  model 
5200  14  channel  tape  recorder.   Real  time  readout  on  an  eight 
channel  chart  recorder,  Brush  model  240,  was  used  to  check 
the  quality  of  the  signals  coming  from  the  sensors. 
4 .   Fluctuation  Statistics 

Spectral  analyses  to  obtain  values  were  performed  on 
fluctuating  velocity  data  obtained  with  Hot  Film  sensors 
oriented  in  the  vertical,  normal  to  the  mean  wind.   These 
data  represent  a  point  measurement  and  yield  temporal  or 
time  descriptions  of  the  fluctuations.   Ten  minute  segments 
of  data  originally  recorded  on  magnetic  tape  were  recorded 
into  an  EMR-Schlumber  model  1510  digital  spectrum  analyzer. 
Procedures  for  converting  spectral  values,  obtained  with  the 
analog  spectral  analyzer,  to  engineering  units  and  for  obtain- 
ing turbulence  parameters  from  the  spectra  are  described  in 
the  following  paragraphs : 

a.   Scaling  Spectral  Plots 

A  necessary  procedure  was  scaling  the  spectral 
plots  to  relate  RMS  input  voltages  to  power  spectral  densi- 
ties (PSD);  variance  per  unit  frequency.   To  obtain  power 
spectral  density  levels,  corresponding  to  RMS  voltage  inputs, 
calibrated  scale  charts  had  to  be  constructed. 

For  purposes  of  the  x-y  plot  format  of  the 
analyzer  output,  the  RMS  voltages  were  converted  to  y  =  log, n 
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(voltage)  units  and  a  graduated  scale  was  constructed  so 
that  the  logarithm  of  volts  RMS  could  be  interpolated  from 
spectral  plots.   The  value  of  the  vertical  scale  (y)  was 
adjusted  for  each  spectrum  as  a  function  of  both  input  gain 
and  spectral  gain.   These  values  were  then  converted  to  PSD 
levels  for  use  in  calculating  e  values. 

The  relation  between  Volts  RMS  to  PSD  units  used 
was 

2 
a/£.    -.c^    .,     (cal.  level  VRMS )     nny>2     /oo> 
S(f)  PSD  units  =  lm5    Bandwidth x  (10*)      (32) 


u  t>   j  fjxv    Freq.  Range 

where       Bandwidth  =  tt — -Vrr i — 

No.  Channels 


,  _  256  Hz    ,  „  » 
(  -  -Ei—  =  1  Hz) 

and  cal.  level  VRMS   =  voltage  at  y  =  0 

(  =  1  VRMS  for  3.16  V  input  setting). 

Amplitude  scaling  calibrations  were  accomplished 
using  externally  generated  "white  noise"  signals  of  1  Volt 
RMS  with  a  frequency  range  from  0.1  Hz  to  10  0  0  Hz  (giving  a 
PSD  of  10"3V2/Hz).   Setting  a  3.16  V  (lOdB)  input  on  the  EMR 
1510  Digital  Spectrum  analyzer  insures  that  1  Volt  corre- 
sponds to  y  =  0.   An  example  of  such  a  calibration  plot  for 
a  frequency  range  of  1-200  Hz  (the  frequency  range  of  inter- 
est) is  shown  in  Figure  10. 

b.   Computation  of  Turbulence  Parameters 
from  Scaled  Spectra 

Values  of  the  turbulence  parameter  e  were  obtained 

from  the  velocity  variance  spectra  on  the  basis  of  the  formula 
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for  the  inertial  subrange  in  wave  number  space  S(k),  Equa- 
tion (6),  which  predicts  a  -5/3  slope  for  the  spectra  when 
plotted  in  a  log-log  format.   Figure  11  is  a  typical  spectra 
considered  in  the  analyses.   Velocity  variance  spectra  con- 
sistently exhibited  slopes  of  -5/3. 

Assuming  -5/3  slopes  for  the  variance  spectra, 
the  intercept  of  the  best  fit  -5/3  slope  with  the  1  Hz  fre- 
quency line  was  the  spectral  density  denoted  (PSD)  value  used 
in  computing  the  parameter  of  interest,  e.   The  measured  PSD 
value  was  converted  to  engineering  spectral  density  units  by 
the  relation 

S  (f)  =  C  2  •  PSD 
u       u 

=  (cm/s/Volt)2-Volt2/Hz  =(cm/s)2/Hz   (33) 


where  C   was  the  calibration  factor  for  the  constant  tempera- 
ture, anemometer.   Determination  of  the  value  for  C   has 

u 

been  described  by  Lund  (1975). 

Since  the  velocity  fluctuations  were  measured  at 
a  fixed  point  in  the  flow,  the  resultant  spectral  values  are 
defined  at  "temporal"  frequencies  denoted  as  S  (f)  in  Equa- 
tion (33). 

5 .   Analysis  of  Paired  Temperature  Sensor 

2 
Cm   values  were  estimated  by  using  paired  platinum 

wire  sensors  for  variance  analysis.   This  analysis  was  based 
on  the  expression  for  the  temperature  structure  function, 
Equation  (5).   The  analysis  is  one  in  which  temperature  differ- 
ences between  the  paired  resistance  wires  separated  by  a 
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distance  r  is  subjected  to  variance  analyses  as  indicated  by 
the  term  (Tf (x)  -  Tf(x+r))2. 

Voltages  corresponding  to  the  temperature  differences, 
AT,  as  measured  by  the  paired  sensors  were  recorded  on  analog 

magnetic  tape.   The  variance  or  root  mean  square  of  the  tem- 

2 
perature  differences  required  to  compute  CT   was  obtained  by 

analyzing  the  analog  voltages  with  a  mean  square  voltmeter 
(Thermo  Systems  Incorporated,  model  1060)  using  a  3  0  second 
time  constant  for  the  averaging.   The  output  of  the  mean-square 
voltmeter  (as  well  as  a  record  of  the  original  signal)  were 
recorded  on  standard  strip  charts. 

Ten  minute  averaged  values  of  the  RMS  voltages  were 
abstracted  from  the  strip  chart  and  these  values  were  utilized 
in  a  program  for  the  IBM  360  computer.   This  program  con- 
verted the  mean  square  voltages  values  to  temperature  differ- 

2 
ence  variance  values,  from  which  C„   were  computed. 

6 .   Fluctuating  Sensor  Calibration 

Platinum  wire  resistance  is  related  to  temperature  by 


R  =  R  (1  +  a(T-T  ) )  (34) 

o  o 

where  R  =  resistance,  T  =  calibration  temperature,  T   =  0  C, 

R   =  resistance  at  T  ,  and  a  =  .0036/  C  (temperature  coeffi- 
o  o '  r 

cient  of  resistance  for  platinum).   Each  wire's  resistance 
was  measured  in  the  laboratory  at  a  temperature  T,  to  deter- 
mine R  . 
o 
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Temperature  fluctuations  can  be  related  to  wire  re- 
sistance fluctuations  by  differentiating  Equation  (3M-) 


or 


dR  =  aR  dT 
o 


RT  =  aR  TT  (35) 


Recorded  voltage  fluctuations  are  then  related  to  re- 
sistance changes  by  the  transfer  function,  H,  in  the  expres- 
sion 

V  =  HR'  (36) 

where  H  was  determined  for  the  bridges  by  varying  the  resis- 
tance of  a  TSI  1056  Variable  Decade  Probe  Resistance  and 
recording  the  corresponding  voltage  changes.   Voltage  vs 
resistance  is  plotted  and  the  slope  of  This  plot  determines 
the  transfer  function,  H.   For  this  experiment  H  =  10  volts/ 
ohms  . 

Finally,  to  relate  recorded  voltage  fluctuations  to 
temperature  fluctuations,  Equation  (36)  is  substituted  into 
Equation  (35)  giving 


T'  =  (1/H  RQ)  V  (37) 


or 


T1  =  CHV  (38) 


where  Ctj  is  the  temperature  calibration  factor  in   C/volt. 
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Velocity  sensors  calibrations  were  accomplished  in  a 
wind  tunnel  using  the  Hot  Wire  Calibrator  TSI  112  5.   In  this 
process,  a  known  wind  velocity  is  blown  across  the  wire  and 
the  corresponding  voltage  output  is  recorded.   The  sensor  wind 
speed  is  given  by 


2 


.x, 


V"  =  aU2  +  b  (39) 


where  V  is  the  hot-film  voltage  output  and  0  is  the  mean  wind 
speed  for  any  given  level.   The  constants,  a  and  b,  are  the 
calibration  curve  slope  and  intercept  respectively.   The  con- 
stants, a  and  b,  were  derived  in  the  lab. 

To  convert  the  voltage  PSD  levels  to  velocity  units 
requires  a  calibration  factor  given  by 

U1  =  CV?  (40) 

where  C  is  the  calibration  factor  in  m/s/volt,  V  is  the 
voltage  fluctuation,  and  U'  the  velocity  fluctuation. 
Differentiating  Equation  (39)  yields 

UT  =  UvtfVa)  V»  (41) 

which  when  substituted  into  Equation  (40)  yields 

C   =  4VU'2/a  (42) 

The  slope  of  the  plotted  curve  yields  the  value  for  a  in  the 
above  equation. 
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C.   ANALYSIS  PROCEDURES 

The  data  were  edited  for  gross  errors  or  inconsistencies 
due  to  known  instrumentation  malfunctions  or  limitations  of 
instruments.   An  example  of  limitation  on  instruments  is  the 
relative  wind  direction  being  greater  than  30°  off  the  bow  of 
the  ship.   Earlier  experiments  have  shown  that  when  these  con- 
ditions exist,  the  winds  are  greatly  influenced  by  the  ship 
and  therefore  are  of  little  value  to  the  experiment. 

The  criteria  for  evaluating  or  discarding  data  from  in- 
dividual levels  or  for  entire  interval  depended  on  performance 
check  results  obtained  during  measurements,  obvious  inconsis- 
tencies, and  the  availability  of  fluctuation  data.   Since  mean 
wind  0,  mean  temperature  T,  and  mean  humidity  q,  are  expected 
to  vary  logarithmically  with  height,  these  parameters  were 
plotted  for  each  time  period  using  the  IBM  360  and  NFS  FORTRAN 
program  PLOTP.   A  best  fit  straight  line  was  drawn  to  the  data 
points.   (See  Figure  12.)   This  procedure  was  very  subjective, 
but  was  intended  to  further  refine  the  data  by  detecting  more 
obvious  discrepancies  and  determining  if  the  data  did  indeed 
vary  logarithmically  with  height.   Gradient  values  of  the 
parameters  used  in  actual  computations  were  evaluated  using  a 
least  square  fit.   By  means  of  the  subjective  procedure,  it 
was  determined  that  the  humidity  sensor  at  the  second  level 
was  reading  too  low  during  most  of  the  experiment.  (See  Figure 
11a.)   This  error  was  corrected  for  the  least  square  fit  by 
assigning  the  second  level  the  same  value  as  level  three. 
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Fig.  12.   Typical  profiles  for  a)  specific  humidity  -  25  Sept 
75  (0105);  b)  temperature  -  25  Sept  75  (0250);  and 
c)  winds  -  25  Sept  75  (0551), 
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IV.   RESULTS 

A.   CT2  RESULTS 

1.   Gradient  Description 

2 
A  comparison  of  36  8  overwater  CT   and  profile  results 

with  the  overland  empirical  expression,  Equation  (22),  appears 
in  Figure  13.   The  distribution  predicted  by  Equation  (22)  is 
the  solid  line.   The  circles  are  mean  dimensionless  tempera- 
ture structure  function  parameter  (DTSFP)  results  over  .50 
Richardson  number  (Ri)  intervals.   The  Standard  deviation 
corresponding  to  each  mean  value  appears  as  vertical  lines 
on  which  the  number  of  values  defining  the  mean  appears.   Sig- 
nificantly, these  results  represent  more  than  a  factor  of  3 
(368  versus  119)  of  those  presented  previously.   (Davidson, 
et  al.,  1976).   However,  the  mean  DTSFP  distribution  changed 
very  little  between  the  two  sets.   Larger  DTSFP  values  for 
the  overwater  regime  for  near  neutral  and  stable  conditions 
(Ri  >  -.5)  can  be  qualitatively  related  to  the  wave  influence. 
Decreased  DTSFP  values  for  unstable  conditions  (Ri  <  -.5)  are, 
perhaps ,  associated  with  aerodynamic  smooth  properties  for 
the  sea  surface  during  light  wind  conditions . 

A  curve  describing  the  observed  data  was  obtained 
using  a  best  fit  polynomial  analysis.   The  results  of  this 
analysis  were 
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=    2.37(.74+Ri)     5#04       Ri    >    0  (43) 


r    2/74/3,86,2 


=    0.56(.24+Ri)~2,13    -1.25    <    Ri    <    0 


=    0.2 


(44) 
Ri    <    -1.25       (45) 
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Fig.  13.   Overwater  results  for  DTSFP  versus  Richardson 
number. 


The  influence  of  the  waves  can  be  further  examined  on 

2     ... 
the  basis  of  observed  height  variations  of  CT  .   Similarity 

2  .        ... 

theory  predicts  that  CT   has  a  height  variation  which  depends 

on  stability  which  was  discussed  by  Wyngaard  et  al.  in  con- 

2  . 
junction  with  an  examination  of  Equation  (12).   Since  T...   is 

independent  of  height  in  the  surface  layer,  Equation  (12) 

yields  the  following  general  predictions  for  the  variation  of 

2 
C„   with  height  based  en  stability. 

2 
a.   For  near  neutral  conditions,  (Z/L  =  0),  CT   de- 

-2/3 

creases  with  height  as  Z      Equation  (12)  becomes 


CT2  =  T,2  Z"2/3 


(46) 
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b.  For  unstable  conditions  (Z/L  <_  0,  Equation  (13)), 

2  -2/3 

CT   decreases  more  rapidly  with  height  than  Z 

2 
If  -7Z/L  >>  1,  Crn      would  decrease  with  height  as 

-4/3 

Z     ,  since  Equation  (12)  becomes 

2/3 

2    „  T  .-    4/'3       -U/3 

CL,      ~    4/3    ^ (wTO    7~    Z    4/ci                    (47) 

T  kg 

2/3  i+/  3 

with  T    /kg(wT')     independent  of  height. 

2 

c.  For  stable  conditions  (Z/L  >  0),Cm   will  decrease 

—      i 

-2/3 

with  height  less  rapidly  than  Z      since 

CT2  =  T:,.  4.9Z~2/3  [1+2.75  h  (48) 


and  the  term  in  brackets  increases  with  increas- 
ing height. 

The  wave  influence  can  also  be  seen  in  height  versus 

2       2  2 

CT   or  CN   (computed  from  CT   values)  results  in  Figures  14 

and  15  respectively.   In  Figure  14,  the  feature  which  suggests 

-2/3 
the  wave  influence  is  the  departure  from  the  Z      height  dis- 
tribution as  the  surface  is  approached.   In  the  upper  three 

2  .      -2/3 

levels  CT   appears  to  fit  a  Z      distribution,  but  as  the 

surface  is  approached  the  distribution  appears  to  be  better 

-4/3 
defined   by  a  Z      height  distribution.   In  Figure  15,  a 

similar  change  is  observed.   The  upper  levels  appear  to  fit  a 

-2/3  2 

Z      distribution  for  C^   but  as  the  surface  is  approached 

-4/3 
fit  a  Z      height  distribution. 

2 
Further  evidence  of  wave  influence  on  CT   appears  in 

joint  probability-conditional  mean  results  of  the  normalized 
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Fig.  14.   Normalized  CT   results  versus  height  (Z),  114 
periods . 


100 


z 

(m) 


10 


10 


a -9/76 
•  -2/77 


C2~Z-4/3 


^N.z^N.ee 


Fig.  15.   Normalized  C„   (computed  from  CT2) 
versus  height.  T 


results 
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temperature  variance  (cJrp/T... )  over  waves  in  Figure  16.   The 

results  and  the  measurements  yielding  them  have  been  de- 

2     2 
scribed  by  Davidson  (1974).   The  ratio  (a^  /T...)   is  direct- 
ly related  to  the  DTSFP,  defined  by  Wyngaard,  et  ai.  (1971). 
In  Figure  16,  the  vertical  axis  variable  (C/U...)  defines  wave 
influence  on  aT/T... ,  C  is  the  phase  speed  of  prevailing  swell, 
L  is  the  Monin-Obukhov  stability  length.   The  results  were 
obtained  under  stable  conditions.   The  diagonal  orientation 
of  the  contours  of  the  number  joint  occurrences  provides 
qualitative  evidence  that  the  waves'  influence  (C/U..c)  on 
<7m/Ta,  hence  the  DTSFP,  is  as  significant  as  the  stability 
influence  (Z/L) . 
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Fig.  16.   Joint  probability-conditional  mean  results  of 
aT/T...  for  C/U...  and  Z/L  dependence. 
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2 .   Bulk  Aerodynamic 

Results  obtained  from  analyses  with  bulk  aerodynamic 
variables,  Equation  (31),  are  compared  to  the  Wyngaard  et  al. 
prediction  curve  in  Figure  17.   The  latter  is  the  solid  line. 
The  circles  are  mean  dimensionless  structure  function  para- 
meter results  over  1.00  X  intervals.   Standard  deviations 
corresponding  to  mean  values  appear  as  vertical  lines  and  the 
number  above  these  vertical  lines  indicates  the  number  of 
values  defining  each  mean.   Larger  values  of  DTSFP  are  again 
attributed  to  the  wave  influence  as  explained  earlier. 

The  following  was  evaluated  for  the  bulk  aerodynamic 
expressions  using  a  best  fit  polynomial  analysis. 


2  2/3      2 

cTzzz/d/(ATr 


=  .037(.642+X)  1,?2   Stable       (49) 


=  .042(.323-X)  *57    Unstable     (50) 


Three  points  from  an  analysis  by  Friehe  (1976)  appear 

in  Figure  17.   A  significant  aspect  of  Friehe 's  analysis  of 

2 
available  data  was  the  exclusion  of  CT   data  from  temperature 

fluctuation  records  which  contained  large  numbers  of  spikes 

exhibiting  ramplike  signatures,  as  shown  in  Figure  18.   This 

led  to  a  significant  reduction  in  the  amount  of  data  used  in 

his  comparison.   Based  on  observations  by  Hicks  (1972)  and 

hypothesis  of  Schmitt  et  al .  (1976),  Friehe  attributed  the 

occurrence  of  the  spikes  to  salt  nuclei  accumulations  on  the 

wires  with  subsequential  erroneously  sensed  temperature 

changes  by  the  system  as  water  vapor  is  absorbed  or  is 

evaporated  on  the  salt  nuclei. 
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Fig.  18.   Time  series  showing  "cold  spikes"  in  the 
temperature  field  under  stable  conditions 
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The  occurrence  and  nature  of  the  spikes  were  also 
examined  with  data  collected  over  Lake  Michigan  by  Davidson 
and  Safley  (1972).   Figure  19  includes  ramp  signatures  in 
temperature  traces  observed  over  Lake  Michigan  where  salt 
accumulation  problem  would  not  be  a  factor.   Therefore,  the 
ramp  signatures  is  possibly  due  to  factors  associated  with 
the  nature  of  the  shear  flow  itself.   On  this  basis  there 
did  not  appear  to  be  sufficient  evidence  to  discard  data 
with  spikes  with  ramp  signatures.   It  is  not  believed  that 
this  data  being  included  in  this  study  caused  the  signifi- 
cant difference  between  the  predicted  Wyngaard  curve  and  the 
observed  data  curve. 

A  contradiction  occurs  between  results  in  Figures  13 

3  6 

and  17.   In  Figure  13,  given  the  same  Z  and  -k—  values  for 

both  regimes,  the  plotted  curves  predict  a  higher  value  of 

2  ... 

C™   for  the  overwater  regime  for  Ri  >  -.5  and  a  higher  value 

2 

of  CT   in  the  overland  regime  for  Ri  <  - . 5 .   On  the  other 

hand  in  Figure  17 ,  given  the  same  Z  and  AT  values  for  both 

2  . 

regimes,  the  plotted  curves  predict  a  higher  value  of  Cm   in 

2 

the  overwater  regime  ror  X  <  1.5  and  a  higher  value  or  C_ 

in  the  overland  regime  for  X  >  1.5.  Since  Ri  and  X  are  cor- 
related they  should  predict  basically  the  same  result.  How- 
ever, as  is  seen  in  Figures  13  and  17,  the  only  ranges  where 
they  correlate  is  the  near  neutral  range. 

In  order  to  explain  the  reason  for  this  contradic- 
tion, values  of  Ri  were  plotted  versus  values  of  X,  Figure  2C 
The  circles  show  average  values  of  Ri  for  X  intervals  of  1  •  • 
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Fig.  19.   Temperature  trace  and  velocity  fluctuation  vectors 
and  wave  heights  associated  with  microthermals 
observed  over  Lake  Michigan. 
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Since  Z/L  is  correlated  with  X,  a  similar  relationship 
should  be  seen  between  X  and  Ri  as  exists  between  Z/L  and  Ri . 
The  depicted  curve  is  similar  to  the  latter  relationship. 
The  conclusion  is  that  the  expected  correlation  existed  be- 
tween these  two  stability  parameters. 

A  second  possible  reason  for  the  contradiction  is  the 

relationship  between  Z    /AT2and  1/Z    (89/3z)2  which  norma- 

2  . 

lize  CT   in  the  DTSFP.   Mean  values  of  the  former  versus  the 

latter  appear  in  Figure  21.   There  should  be  a  positive  cor- 
relation between  the  two  and  this  occurs  only  for  large  values 

4/3        2 
of  1/Z    (99/8z)  .   For  small  values  of  the  latter  there  is, 

however,  a  negative  correlation.   This  result  indicates  thar 
the  two  parameters  are  correctly  correlated  when  small  tem- 
perature gradients  exist,  but  are  not  when  the  gradients  be- 
come large.   The  contradiction  is  likely  to  be  due  to  the 
•problem  of  accurately  defining  the  bulk  temperature  difference, 
AT. 

B.   s  RESULTS 

1.   Gradient  Description 

Comparison  of  the  observed  £  results  and  the  curve 
based  on  Equations  (23)  and  (24)  is  shown  in  Figure  2  2  where 
the  predicted  curve  is  the  solid  line.   Averages  over  Ri  in- 
tervals of  .25  of  observed  data  appear  as  small  circles. 
The  error  bars  are  standard  deviations  from  the  mean  within 
each  interval.   The  number  of  observations  defining  each 
mean  value  is  at  the  top  of  the  error  bars.   The  significant 
difference  between  the  predicted  curve  and  the  observed  data 
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for  Ri  >  .20  is  expected  on  the  basis  of  Equation  (23).   For 
Ri  >  .20  the  equation  is  no  longer  valid  since  terms  in  the 
expression  become  infinite. 

An  examination  of  Equations  (13),  (14),  and  (15) 

3 
with  the  knowledge  that  U...   and  L  are  independent  of  height 

in  the  surface  layer  yields  the  following  general  prediction 

for  the  variation  of  e   with  height: 

a.  For  neutral  conditions,  (Z/L  =  0),  £  will  decrease 
with  height  as  Z~   since  Equation  (15)  becomes 

£  =  -rr-   Z  x  (51) 

b.  For  stable  conditions,  (Z/L  >  0 ) ,  £  will  decrease 
less  rapidly  with  height  than  Z    since  Equation 
(15)  becomes 

s  =  -— ■  Z  "L(l+3.7  Z/L)  (52) 

c.  For  unstable  conditions,  (Z/L  <  0),  £  will  de- 
crease  less  rapidly  with  height  than  Z  .  since 
Equation  (15)  becomes 

u,3 

Z  ^(1-15  Z/L)~*  -Z/L      (53) 


K 

The  quantity  in  brackets  increases  with  increas- 
ing height.   This  increase  is  less  than  for 
stable  conditions.   So  the  decrease  of  £  with 
height  for  unstable  conditions  is  less  than  that 
for  stable  conditions. 
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Results  in  Figures  23  and  24  indicate  that  under  near  neutral 
conditions  e  fits  a  Z~   height  distribution  quite  well. 

There  appears  to  be  no  change  in  the  height  distribution  as 

o 

the  surface  is  approached  as  there  was  with  Cm  . 

2 .   Bulk  Aerodynamic 

Results  obtained  from  analyses  with  bulk  aerodynamic 
variables  appear  in  Figure  25.   The  predicted  curve  is  the 
solid  line,  while  the  circles  are  mean  non-dimensional  e  re- 
sults over  X  intervals  of  1.0.   Again,  standard  deviations 
corresponding  to  mean  values  are  vertical  lines  and  the 
number  above  these  lines  indicate  the  number  of  values  defin- 
ing each  mean.   The  correlation  between  the  non-dimensional 
s  values  versus  X,  as  described  for  Z/L  by  Busch  and  Panofsky 
(1969),  is  quite  good. 

C.   C  2  RESULTS  AT  REGIONS  ABOVE  SURFACE  LAYER 

A  Cessna  was  used  in  this  experiment  to  extend  the  data 

base  beyond  the  height  of  the  instrumented  towers.   Height 

2 
distribution  of  CT   values  as  expressed  m  Equations  (46), 

(4-7),  and  (45),  have  been  shown  to  be  valid  to  height  above 

-2/3 
4000  meters  (Tsvang,  1969).   Figure  26  shows  that  a  Z 

distribution  exists,  indicating  stable  or  near  neutral  condi- 
tions, from  lowest  level  to  about  100  meters.   Above  100 

2 

meters,  there  is  a  rapid  increase  of  C™   values. 

Frisch  and  Ochs  (1975)  suggested  that  if  the  height  of 
the  inversion  base,  Zi,  is  another  relevant  scaling  parameter, 
then  proceeding  in  a  manner  similar  to  that  of  Wyngaard  et  al, 
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Fig.  26.   CT   height  distribution  and  temperature  profile 
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Equation  (47)  could  be  modified  with  an  additional  non-dimen- 
sional function,  GCZ/Z-).   Equation  (47)  then  becomes 

CT2  =  [_^7T](T/g)2/3(w-TTr)4/3Z-l+/3  S(|  )   (54) 
3k  i 

where  G(Z/Z.)  =  1  for  Z/Z.  <<  1  in  order  to  have  Equation 

(54)  reduce  to  Equation  (47)  in  the  surface  layer. 

-4/3 
Frisch  and  Ochs  found  a  significant  departure  from  a  Z 

law  for  values  of  Z/Z-  >  0.1.   This  indicates  that  the  atmos- 
phere is  no  longer  in  the  "free  convection"  regime,  but  is 
affected  by  the  lid  on  the  convective  layer  at  Z . .   This 
effect  would  restrict  the  integral  scale  of  the  turbulent 
temperature  fluctuations  with  height  compared  to  the  free  con- 
vection value.   In  addition,  because  of  the  layer  scale 
circulation  and  mixing,  the  temperature  variance  should  be 

greater  than  the  free  convection  variance.   These  effects 

2 
will  increase  the  value  of  CT   relative  to  the  free  convec- 

2 
tion  prediction.   Figure  26  shows  the  rapid  increase  of  CT 

values  as  the  inversion  is  approached. 
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V.   CONCLUSIONS 

Application  of  formulated  overland  expressions  in  the 

overwater  regime  showed  a  significant  difference  between  ob- 

2 

served  CT   overwater  data  and  overland  predictions.   This 

difference  is  attributed  to  the  wave  influence  for  near 
neutral  and  stable  conditions  and  to  the  aerodynamic  smooth 
properties  of  the  sea  surface  during  light  wind  periods  for 
moderate  unstable  conditions.   The  formulation  of  these  over- 
water  curves  will  greatly  enhance  the  Navy's  operational 
ability  to  describe  the  marine  boundary  layer  and  its  in- 
fluence on  the  optical  wave  propagation. 

In  comparing  the  DTSFP  overland  and  overwater  when  plotted 
against  Ri  and  then  versus  X,  significant  differences  were 
noted.   Conclusive  explanations  for  this  difference  cannot 
be  given  in  this  study.   Therefore,  recommend  that  serious 
consideration  be  given  to  an  in-depth  examination  of  this 
significant  feature. 

Correlation  of  non-dimensional  e  values  with  the  stability 
parameter,  Ri ,  showed  agreement  between  results  and  predic- 
tions.  For  Ri  >  .20,  the  theory  is  not  valid  and  therefore 
there  was  very  poor  correlation.   Any  attempt  to  describe 
the  normalized  e  for  Ri  >  .20  in  the  future  will  require  an 

expression  which  will  remain  valid  for  all  values  of  Ri . 

2 
CT   results  appear  to  be  described  by  the  same  high  dis- 
tributions up  to  levels  well  above  the  surface  layer, 
approximately  60%  of  the  distance  to  the  marine  inversion. 
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